Common Envelope Evolution 



Ronald E. Taam, ^ Paul M. Ricker ^ 



O 

o 

(N 
> 

o 



> 

o 



o 

o 
U 



X 



^Department of Physics and Astronomy, Northwestern University, Evanston, IL 

60208 

^Department of Astronomy, University of Illinois, Urbana, IL 61801 



Abstract 

The common envelope phase of binary star evolution plays a central role in many 
evolutionary pathways leading to the formation of compact objects in short period 
systems. Using three dimensional hydro dynamical computations, we review the ma- 
jor features of this evolutionary phase, focusing on the conditions that lead to the 
successful ejection of the envelope and, hence, survival of the system as a post com- 
mon envelope binary. Future hydrodynamical calculations at high spatial resolution 
are required to delineate the regime in parameter space for which systems survive 
as compact binary systems from those for which the two components of the sys- 
tem merge into a single rapidly rotating star. Recent algorithmic developments will 
facilitate the attainment of this goal. 
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1 Introduction 



The evolutionary study of close bi- 
nary systems with compact neutron 
star and black hole components has 
been a major focus of stellar X-ray 
astronomy ever since the seminal 
contributions by van den Heuvel and 
Heisedli). Since then, the discovery of 
compact stars in interacting binary 
systems has had an enormous impact 
on our basic understanding of the 
properties of these stars, on the man- 
ner in which mass is transferred from 
a stellar donor to its compact accre- 
tor, and on the evolution of these 



systems (jij; 0). The recent Chandra 
and XMM-Newton observations of 
the point X-ray source population in 
the Milky Way and in external galax- 
ies have led to a renewed interest in 
studying the evolution of these com- 
pact systems in a variety of galactic 
environments. 



Among the trademarks of the evo- 
lutionary scenarios presented in van 
den Heuvel and Heise (jl|) are the car- 
toons used to depict the evolution 
of binary systems to the compact 
stage. These have now become the 
standard means of visualizing the 
possible evolutionary channels. For 
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a recent review of the formation and 
evolution of compact X-ray sources, 
see Tauris & van den Heuvel (jsl). 
Central to the construction of these 
channels is the existence of a com- 
mon envelope or spiral in phase in 
which significant mass and orbital 
angular momentum are lost from the 
system. In this phase, the system 
is transformed from one of long or- 
bital period to one of short orbital 
period (0). One of the members of 
the system, a star with an evolved 
core in the red giant or asymptotic 
giant branch star, expands beyond 
the orbit of its companion. Provided 
that the stellar components are not 
in state of synchronous rotation at 
the onset of mass transfer, the asyn- 
chronous interaction of the two com- 
ponents can lead to the conversion of 
orbital energy into the kinetic energy 
of outflow of the common envelope. 
As a result, it was envisioned that a 
main sequence-like companion spi- 
rals toward the core of the giant. The 
common envelope is ejected leav- 
ing behind the progenitor remnant 
of the compact component with its 
companion at separations which are 
significantly smaller than the orig- 
inal radius of the giant progenitor 
(see reviews by Iben & Livio (0) 
for evolutions involving intermediate 
mass stars and Taam & Sandquist 
(0) for evolutions involving massive 
stars). Among the classes of binary 
systems for which the common enve- 
lope phase plays an important role 
are systems with neutron star (bi- 
nary radio pulsars. X-ray binaries), 
black hole (X-ray sources), or white 
dwarf (cataclysmic variables) com- 
pact components. 



In the next section, the evolution 
leading to the common envelope 
phase is described without reference 
to the evolution resulting from the 
collisional interactions in dense stel- 
lar systems (0; S [o]). Of particular 
importance, here, are the conditions 
leading to the establishment of a 
non-corotating common envelope in 
primordial binary systems. The var- 
ious stages of the ensuing spiral in 
phase are presented in section 3, and 
the conditions leading to the success- 
ful ejection of the common envelope 
are described in section 4. Finally, 
we provide an overview of the results 
from hydrodynamical simulations 
and discuss the prospects of high 
spatial resolution calculations in the 
near future in the last section. 



2 Evolution to the common en- 
velope stage 



The binary system can enter into 
the common envelope stage when its 
total systemic angular momentum 
is less than its minimum value for 
components in synchronous motion. 
In this case, it evolves away from 
synchronism as a result of a tidal 
instability (0; [ill; H H Hi). This 
critical limit for synchronous sys- 
tems exists because the spin angular 
momentum of the more evolved and 
significantly larger stellar compo- 
nent can become comparable to the 
orbital angular momentum of the 
system. A synchronously rotating 
star which evolves to this point can- 
not evolve to ever larger radii and 
still remain synchronized to the or- 
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bital motion. Instead, the two com- 
ponents of the system are driven 
out of synchronism by the action of 
tidal forces and viscous dissipation, 
shrinking the orbit as orbital angu- 
lar momentum is converted into spin 
angular momentum of the larger 
component of the system. For a star 
of a given radius, the spin angular 
momentum {Jspin « a~'^/^) increases 
with decreasing orbital separation, 
a, while the orbital angular momen- 
tum decreases {Jorb on a^^^). Hence, 
for sufficiently small orbital separa- 
tions, the rotational angular velocity 
of the more evolved stellar compo- 
nent cannot be maintained at the 
orbital angular velocity of its com- 
panion. This evolution leading to a 
spiral in phase takes place for binary 
systems characterized by large mass 
ratios. For example, systems with a 
red giant or asymptotic giant branch 
component would undergo this tidal 
instability for mass ratios greater 
than about 5-6. 

An alternative path to the common 
envelope stage can occur whenever 
the rate of mass transfer is suffi- 
ciently high such that the timescale 
for angular momentum redistribu- 
tion between the spin and the orbit 
required to bring the components 
of the system into a state of syn- 
chronous rotation is longer than the 
mass transfer timescale. This evolu- 
tionary channel can take place when 
the mass transfer rate accelerates 
due to an imbalance between the 
variation of the stellar radius and the 
star's corresponding Roche lobe with 
respect to mass loss. This mass trans- 
fer instability is particularly relevant 
to systems in which the more massive 
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Fig. 1. The orbital separation as a func- 
tion of time for a system composed 
of an asymptotic giant branch star of 
3Mq with a core mass of 0.7Mq syn- 
chronously rotating with a OAMq main 
sequence companion for a system with 
an initial orbital period of 0.84 yr (28). 



star is characterized by a deep con- 
vective envelope. In this case, there 
is a tendency for the donor to expand 
while its corresponding Roche radius 
decreases as mass is lost. The accre- 
tor is also likely to expand since it 
accretes at a rate faster than its ther- 
mal timescale and cannot assimilate 
the accreted mass in a state of ther- 
mal equilibrium. As a result, both 
components are expected to fih or 
exceed their respective Roche lobes. 
In this pathway, systems in which the 
more massive star undergoes mass 
transfer while in the Hertzprung gap 
or in the red giant or asymptotic 
branch phase can enter the common 
envelope phase either via a direct 
or a delayed dynamical instability 
(jlSl ). We note that evolution to the 
common envelope stage via a mass 
transfer instability is not limited to 
accreting non degenerate stars, but 
can also apply to compact accreting 
neutron stars ( 1l6l ) as well. 
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3 Evolutionary stages 



Substantial progress in the descrip- 
tion of the spiral in evolution has 
been achieved from investigations 
of increasing complexity and real- 
ism. Building ^on the early studies 

19l ) and two 



we now turn. 
1 



in both one 
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( 1201 : l2ll : l22l : |23| ) dimensional model- 
ing, the studies based on three di- 



mensions (124; l27|; l25|; M, l28|; l29|) have 
provided significant understanding 
of the conditions leading to success- 
ful ejection of the common envelope. 
A particularly challenging aspect of 
the mult i- dimensional calculations 
is the wide range of timescales and 
length scales that must be modelled. 
In particular, the timescales range 
from hours during the late stages to 
years at the onset of the common en- 
velope phase, and the length scales 
range from Earth-like dimensions to 
more than an astronomical unit. As 
a result, the hydrodynamical calcu- 
lations that have been reported are 
primarily exploratory in nature. On 
the other hand, the existing com- 
putations have provided insight into 
the essential features introduced by 
the additional dimensionality con- 
sidered in the problem and the im- 
portant physical processes required 
for a detailed description underlying 
our understanding of this evolution- 
ary channel. In the following, we 
provide an overview of the several 
distinct stages of the common enve- 
lope phase, focusing on the evolution 
during the deep spiral in phase to 
short orbital periods and the termi- 
nal phase leading to the formation 
of the final post common envelope 
system. To these individual phases 
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Fig. 2. The density distributions and ve- 
locity fields in the orbital plane of the 
binary described in Figure [1] (|2^). The 
density contours are logarithmic and 
are spaced five per decade. The velocity 
vectors are scaled to the maximum in 
each panel. The solid dots indicate the 
position of the cores and the distribu- 
tions are at times given in each panel. 
Spin up can be seen at a time of 1.02 
yrs with material outflow subsequently 
taking place (e.g., see panels at times 
1.26 and 1.45 yrs). 



3. 1 Deep spiral in 



To illustrate a typical binary system, 
in the common envelope phase, in 
Figure [1] we show the orbital sepa- 
ration function of time for an 
asymptotic giant branch star with 
a low mass main sequence compan- 
ion. It is clearly evident that the 
separation between the two cores 
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decreases rapidly (on a timescale of 
about 100 days). During the inspi- 
ral, the outer layers of the common 
envelope are tidally stripped in the 
orbital plane, leading to the forma- 
tion of an outflowing circumbinary 
disk of material. The relative orbit of 
the two cores becomes very eccentric 
during the early phase of evolution 
as most of the orbital angular mo- 
mentum is converted to spin angular 
momentum in the outer envelope 
layers. However, as the two cores spi- 
ral closer together, a greater fraction 
of orbital energy is lost at smaller 
separations in comparison to orbital 
angular momentum, causing the or- 
bit to become more circular. 



The rapid orbital decay of the system 
causes envelope material to acceler- 
ate relative to the companion star, 
leading to the formation of shocks in 
the envelope. As the system evolves 
to ever smaller orbital separations, 
the shocks evolve into a tighter and 
tighter spiral structure pattern. The 
interaction of the shocks with the en- 
velope converts the orbital angular 
momentum of the two stellar compo- 
nents into the spin angular momen- 
tum of the common envelope. As a re- 
sult, matter in the common envelope 
is spun up to nearly circular motion 
(see figure[2]). This spin up effectively 
causes an outward force, as centrifu- 
gal effects help to expand material to 
greater distances from the two cores. 
As shown in Figure [3l common en- 
velopes involving massive supergiant 
stars also show this behavior m). 




X (It)'' cm) 



Fig. 3. The density distributions and 
velocity fields in the orbital plane of a 
binary composed of a 20Mq red super- 
giant with a 1.4 Mq companion at dif- 
ferent evolution times (@). Density con- 
tour levels correspond to 5 per decade, 
and the velocity field is scaled to 60 km 
s~^ in all panels. The evolutionary state 
of the massive star is at the onset of 
core carbon burning. A tightly wound 
spiral emerges as the orbital velocities 
of the two cores exceed the rotational 
velocity of the envelope gas. The solid 
dots indicate the position of the cores 
and the initial orbital period is 10.6 yr. 
Spin up occurs at time 10.4 yr, leading 
to an outward flow at a time of 12.7 yr. 

3. 2 Final stage 

Since the orbital energy is deposited 
into the common envelope on a 
dynamical timescale in the orbital 
plane, the ejected material is concen- 
trated toward this plane. The ejec- 
tion, however, is nearly axisymmet- 
ric with respect to the angular mo- 
mentum axis since the timescale on 
which the energy is deposited (com- 
parable to the local orbital period) 
is less than the timescale of orbital 
decay in this stage. The flow is not 
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Fig. 4. The density distributions and ve- 
locity fields in the plane perpendicular 
to the orbital plane at the times for the 
system in Figure [2] (0). Density con- 
tours and velocity vectors are as in Fig- 
ure [21 

purely radial, as can be seen most 
clearly at a time of 1.45 yr in Figure 
m where a circulatory flow between 
the equatorial plane and the polar 
direction results. 

Spin up to a greater degree is found 
for a red giant progenitor with a low 
moment of inertia either due to low 
envelope mass or relatively small 
radius. In this greater den- 

sity contrast between the polar and 
equatorial directions can develop. 
This may be relevant to the bi-polar 
morphology of a class of planetary 
nebulae. Figure [5] illustrates the den- 
sity and velocity distribution in the 
plane perpendicular to the orbital 
plane for two systems in which the 
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Fig. 5. The density contours and veloc- 
ity field in the plane perpendicular to 
the orbital plane after the rapid infall 
phase (jjgl ). Density contours are 5 per 
decade, and the velocity fields in each 
panel are scaled to a maximum speed 
of 60 km s~^. The mass of the red giant 
core and companion in the common en- 
velope binary are O.45M0 and O.SSMq, 
respectively. The mass of the red giant 
in panel (a) is IMq, and in panel (b), 
2Mq. Note that a sharper density con- 
trast between equatorial and polar di- 
rections results in the case for which the 
mass of the companion is comparable to 
the mass in the common envelope. 

mass of the companion and the mass 
of the red giant core are O.SSMq and 
0.45Mq respectively. The structure 
in panel (a) for a IMq red giant and 
in panel (b) for a 2Mq red giant illus- 
trates the dramatic difference in the 
common envelope structure when 
the degree of spin up is great. Note 
the evacuation of matter along the 
polar directions. Eventually an out- 
flow is produced that is concentrated 
in the equatorial plane. 

The evolution enters the terminal 
phase as the companion further spi- 
rals into the common envelope. The 
orbital decay rate decelerates be- 
cause most of the common envelope 
material has moved to larger radii, 
where its interaction with the two 
cores is greatly diminished. As a 
result, the orbital decay timescale 
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Fig. 6. The orbital decay timescale as 
a function of time corresponding to the 
temporal variation of the orbital sepa- 
ration in Figure [1] (0) . 

dramatically increases from about 
100 days to more than 10 years (see 
Figure [6]) . 

At this point we expect that the ma- 
terial interior to the orbit will con- 
tract on its local thermal timescale to 
the degenerate core of the giant. Al- 
though this phase has not been fol- 
lowed, a post spiral in system should 
be formed with the companion star in 
close proximity to the remnant core 
of the former giant. 



4 Conditions for successful 
ejection of the common en- 
velope 

The hydrodynamical simulations 
demonstrate that material in the 
common envelope is ejected at veloc- 
ities greater than that required for 
escape. Although the calculations 
have not been evolved to the phase 
in which the entire envelope is lost 
before the two cores merge, the re- 
sults to date do provide insight into 
the conditions necessary for success- 



ful ejection of the common envelope 
to take place. 

The primary requirement for ejec- 
tion of the common envelope is based 
on energetics. Namely, the energy re- 
leased from the orbit must exceed the 
binding energy of the common enve- 
lope from both stellar components of 
the system. Generally, the efficiency 
of mass ejection, as measured by the 
ratio of the binding energy to the en- 
ergy lost from the orbit, is less than 
unity. 

Secondly, the rate at which energy 
is lost from the orbit must be suffi- 
ciently high that it can be directly 
converted into the energy of the out- 
flow rather than being transported 
to the common envelope surface, 
where it can be radiated away. For 
the dynamical evolution described 
above, the ejection process is rapid 
and, hence, adiabatic. However the 
efficiency of mass ejection is less than 
unity since the mass in the common 
envelope is preferentially ejected in 
the orbital plane of the binary sys- 
tem at velocities which are greater 
than the escape speed. 

Finally, the timescale on which the 
mass is lost from the system must 
be shorter than the inspiral timescale 
of the binary, for otherwise the two 
cores will spiral together and merge 
before the entire common envelope is 
ejected. 

This latter condition points to the 
importance of a core envelope struc- 
ture characteristic of evolved stars 
on the giant branch. For such stars, 
the density profile above the nuclear 
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burning shells is steep so that the 
mass enclosed in an extensive region 
is small. Examples of this particular 
structure are illustrated in Figure [7] 
for a 2OM0 star in different stages of 
core helium burning and in Figure [8] 
for a IMq red giant star character- 
ized by a range of degenerate helium 
core masses. 
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Fig. 7. The variation of mass (in Mq) 
with radius for a 20Mq red super- 
giant. The dashed, dashed dot, and sohd 
curves correspond to a central hehum 
content of 0.5, 0.34, and respectively 
(@). Note the extensive region between 
5 X 10^° and 5 x 10^^ cm in which the 
mass varies very gradually with radius 
at the stage of carbon core burning. 

In both figures an extensive region 
containing little mass is present. The 
fiat mass profile naturally develops 
in the late core helium burning stage 
of a massive star when the central 
helium content nearly vanishes. For 
more advanced evolutionary states, 
the size of this region increases. 
Similarly, red giant stars show such 
profiles for helium degenerate core 
masses greater than about O.2M0, 
extending beyond a solar radius for 
helium core masses greater than 
about O.27M0. 




Fig. 8. The mass profile for a IMq 
red giant star in different phases of 
evolution (fiil ). The solid, dashed-dot, 
and dashed curves denote evolutionary 
phases in which the mass of the degen- 
erate helium core is 0.19, 0.23, and 0.27 
Mq respectively. 

In cases for which the loss of orbital 
energy does not lead directly to hy- 
drodynamical ejection of the deep 
layers of the common envelope, mass 
can still be lost, but on a timescale 
significantly greater than the dynam- 
ical timescale. Here, the rate of en- 
ergy loss from the orbit supplements 
the stellar luminosity generated by 
the nuclear burning shells, enhanc- 
ing the rate of mass loss via the 
radiative processes responsible for 
stellar winds. Such a circumstance 
may lead to the ejection of the com- 
mon envelope, especially when the 
envelope is sufficiently spun up. De- 
tailed calculations reveal that such 
conditions occur when the mass of 
the common envelope is comparable 
to the mass of the inspiraling compo- 
nent. This outcome for the common 
envelope phase may be related to an 
alternative prescription for common 
envelope evolution based on angu- 
lar momentum arguments under the 
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implicit assumption that energy is 
conserved, (jsol ). 

Independent of whether the final 
phase of the common envelope is 
described dynamical or more 

gradual event, the successful ejection 
of the envelope requires progenitor 
stars in the red giant or asymptotic 
giant phase. For progenitors on the 
main sequence or in the Hertzsprung 
gap, it is likely that the binary com- 
ponents will merge in the process, 
leaving behind a single rapidly rotat- 
ing star, with possibly a circumstel- 
lar disk. The class of rap idly rotating 
FK Comae giants (|3l!) may be an 
observational manifestation of such 
a merger event. 



5 Summary and future work 

Computational hydrodynamical in- 
vestigations of the common envelope 
phase have led to a number of re- 
sults with significant bearing on the 
outcome of the post spiral in state of 
the binary system. 

5. 1 Overview of results 

• The binary orbit decays rapidly af- 
ter the onset of the common enve- 
lope phase, occurring on a dynam- 
ical timescale (i.e., on a timescale 
comparable to the orbital period 
for which the system entered the 
common envelope phase). 

• As a result of the effectiveness of 
the gravitational torques in remov- 
ing orbital angular momentum 
from the system, the orbit of the 



binary system becomes eccentric. 
Circularization of the binary orbit 
takes place in the deep layers of 
the common envelope as a signifi- 
cant fraction of the orbital energy 
is lost during this phase. 

• There is strong evidence for signifi- 
cant spin up of gas surrounding the 
region containing the core of the gi- 
ant progenitor and the inspiralling 
companion, reducing the effective 
gravity and making mass ejection 
easier. 

• Matter is ejected from the com- 
mon envelope in all directions, 
with a preference for the orbital 
plane of the binary system, pro- 
ducing a density contrast in the 
ejecta between the equatorial and 
polar directions. 

• The efficiency of the mass ejection 
process, as measured by the ratio of 
the binding energy of the envelope 
with respect to the binary system 
to the orbital energy released from 
the binary during the common en- 
velope phase, is less than about 40- 
50%. 

• The stabilization of the orbit at 
small separations during the final 
phase of the common envelope is 
easiest for progenitor stars which 
have steep density gradients above 
the nuclearly evolved core (so that 
the mass profile is fiat) and/or 
where the mass of the common en- 
velope is comparable to the mass 
of the inspiralling companion. 

Based on these results, one finds that 
the survival of the remnant binary as 
a post spiral in system resulting from 
the successful ejection of the common 
envelope is favored for systems in 
which the giant-like progenitor star 
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is characterized by larger ratios of 
the core mass to total mass (roughly 
greater than about 0.2). Thus, for- 
mation of cataclysmic variable type 
systems is expected to be viable for 
red giant and asymptotic giant star 
progenitors with core masses greater 
than about 0.2Mq and O.GMq re- 
spectively. For more massive stellar 
progenitors (i.e., for stars more mas- 
sive than about 12Mq) relevant to 
the formation of low mass X-ray bi- 
nary and intermediate mass X-ray 
binary systems, ejection of the com- 
mon envelope is favored for advanced 
evolutionary stages of the progenitor 
star during its late core helium burn- 
ing stage and beyond. We note that 
for stars more massive than about 
40Mq stellar winds may significantly 
affect their evolution, precluding the 
Roche lobe overflow and mass trans- 
fer phase since such stars do not 
enter the red supergiant phase (1321 ). 
For systems entering the common 
envelope phase, significant shrink- 
age of the binary orbit is possible, 
leading to a reduction of the orbital 
separation by more than a factor of 
100 from its initial value. 



5.2 Ongoing work 



Although the calculations that have 
been carried out in previous stud- 
ies have provided much insight into 
the phases of the common envelope 
stage, calculations at high spatial 
resolution will be necessary to fur- 
ther quantify the outcome of the 
common envelope phase. In recent 
years the development of sophisti- 
cated computer methodologies has 



made it possible to achieve this goal. 
Specifically, adaptive mesh refine- 
ment techniques have advanced to 
the point where such calculations 
can now be envisioned. For a re- 
cent review see Norman (jssl ). These 
methods allow one to maintain high 
spatial resolution in the core regions 
as the two cores spiral in together. 




Fig. 9. The orbital separation clS Si func- 
tion of time for a system composed of 
a red giant branch star of IMq with 
a 0.7Mq main sequence companion in 
which the initial orbital period is 1 
month. 

Recently, we have initiated high 
resolution calculations using such 
techniques. Since the core of the 
giant and its companion are both 
too compact to resolve even with 
an adaptive mesh, and because of 
their much higher densities with re- 
spect to the matter in the common 
envelope, they couple to the gas pri- 
marily via their gravitational fields. 
Therefore we have modelled the gi- 
ant core and its companion as par- 
ticles using an body solver based 
on the particle mesh method (js^)- 
The companion is represented by a 
single particle, and the giant core 
is represented by a uniform, spher- 
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ical cloud of 2 x 10^ particles with 
core radius three times the smallest 
zone spacing. All of the particles in 
the cloud move rigidly together with 
the cloud's center of mass. This ar- 
rangement ensures that the mapping 
of particle densities onto the mesh 
and of the gravitational forces onto 
the cloud's center of mass are free 
of Cartesian grid effects introduced 
by the use of a cloud-in-cell mapping 
kernel. An adaptive multigrid solver 
is used to solve the Poisson equa- 
tion for the gravitational potential 
(1351: l36l) . Multigrid algorithms accel- 
erate the convergence of relaxation 
methods by covering the computa- 
tional domain with a hierarchy of 
meshes with different spacings and 
applying relaxation to each mesh. 
The single-mesh convergence rate is 
controlled by error modes having the 
longest wavelengths compared with 
the mesh spacing. Thus the longest 
wavelengths in the domain con- 
verge most rapidly on a very coarse 
mesh. By combining results from all 
meshes a multigrid algorithm brings 
all wavelengths into convergence at 
the same rate. 

The use of the adaptive mesh re- 
finement technique and the multi- 
grid method for the gravitational 
potentials improves on our previous 
studies ( l28l : 1291 ). allowing calcula- 
tions not limited to systems of ex- 
treme mass ratio, in which the center 
of mass moves only slightly during 
the evolution. We have performed 
several preliminary calculations of 
the common envelope phase using 
these techniques in the FLASH code 
(ISTl ). FLASH is an adaptive, parallel, 
multi-dimensional hydrodynamical 



simulation code in which the equa- 
tions of compressible gas dynam- 
ics are solved using the Piecewise- 
Parabolic Method fi38i) on a struc- 
tured, rectangular grid. PPM was de- 
veloped to provide very accurate so- 
lutions for flows containing sharp dis- 
continuities, and it greatly improves 
upon the hydrodynamics solver used 
in the nested grid simulations. 

As an example, we illustrate the com- 
mon envelope evolution of a binary 
system consisting of a IMq red gi- 
ant progenitor with a 0.7 Mq main se- 
quence companion at an evolution- 
ary phase when the dynamical evolu- 
tion occurs on a timescale of about 1 
month. The orbital separation of the 
two cores within the common enve- 
lope binary system is illustrated 
function of time in Figure [91 For the 
case in which the red giant was rotat- 
ing at half the synchronous rate, the 
orbit rapidly decayed by a factor of 2 
within about 1 month, with the rela- 
tive orbit passing through two phases 
of apastron and periastron passage. 

In Figure [101 we show the early 
phases of the spiral in process cal- 
culated using FLASH at three time 
slices at the onset of the evolution 
(upper 3 panels) and toward the end 
of the calculation (lower 3 panels) in 
the innermost 5 x 10^^ cm of the bi- 
nary. The evolution is similar to that 
described in Section 3, with the de- 
velopment of spiral shock waves em- 
anating from the two cores. However, 
most importantly, it is seen that the 
adaptive mesh refinement technique 
performs very well in resolving the 
vicinity of the red giant core and the 
outer layers of the common envelope. 
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Fig. 10. The density distribution in the orbital plane during the initial and spiral 
in phase of a binary system composed of a red giant branch star of IMq with a 
0.7 Mq main sequence companion for a system with an initial orbital period of 1 
month. Each of the six panels provides the evolution time, which ranges from 1.55 
days to 41.52 days. For convenience, a scale of 10^^ cm is also indicated. The black 
dot indicates the position of the red giant companion. 



We remark that such an evolution 
could not have been simulated us- 
ing our previous stationary nested 
grid technique due to the hmitations 
on the mass of the companion star 
which, for sufficiently large masses, 
caused the core of the red giant to 
move significantly off the highest 
spatially resolved grid. 

Evolutionary calculations such as 
these in which the innermost regions 
of the common envelope surrounding 
the two cores are highly resolved will 
be necessary to determine whether 
the binary orbit stabilizes or con- 
tinues to shrink. Such investigations 
are essential for determining the 
dependence of the efficiency of the 
mass ejection process on the mass of 



each component, their evolutionary 
stage, and orbital separation. Using 
our model binary systems we plan to 
systematically investigate the mass 
ejection process as a function of these 
system parameters. This is critically 
important since the generality of the 
hydrodynamical results with regard 
to binary system parameters must be 
ascertained. As two possible mech- 
anisms have been identified for the 
emergence of the system from the 
common envelope phase, an explo- 
ration of their different dependencies 
on the binary system parameters is 
required before definitive conclusions 
can be drawn for use in population 
synthesis calculations. This study 
will allow us to distinguish the re- 
gions in progenitor parameter space 
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for which the post common enve- 
lope systems survive from that in 
which the systems merge into a single 
rapidly rotating object. Such results 
will provide important input for pop- 
ulation synthesis modeling, allowing 
us to determine the birth rates for 
the formation of compact systems 
in a variety of galactic environments 
(metal rich star burst regions vs. 
metal poor regions). The comparison 
of such rates with observations will 
allow us to establish the evolutionary 
channels for these celestial sources. 
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